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Abstract 

Using variable temperature scanning tunneling microscopy and dl/dz barrier height spec- 
troscopy, the structure of water on /j-BN/Rh(lll) nanomesh has been investigated. Below 
its desorption temperature, two distinct phases of water self-assemble within the 3.2 nm unit 
cell of the nanomesh. In the 2 nm holes, an ordered phase of nano-ice crystals with about 
40 molecules is found. The ice crystals arrange in a bilayer honeycomb lattice, where the 
hydrogen atoms of the lower layer point to the substrate. The phase on the 1 nm wires, is a 
low density gas phase, which is characterized by contrast modulations and streaky noise in 
the STM images. Tunneling barrier measurements infer the proton positions in the nano-ice 
clusters. 
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Introduction 



Boron Nitride Nanomesh: A template . 



Templates are essential in nanoscience and technology, since they allow the self-assembly of seem- 
ingly artificial structures with new functionalities. If a surface shall act as a natural template that 
imposes a lateral ordering on a length scale larger than a 1 x 1 unit cell, it has to rely on recon- 
struction and the formation of superstructures. For the case of the /z-BN/Rh(lll) nanomesh^ the 
superstructure bases on 12x12 Rh unit cells on top of which 13 x 13 BN unit cells coincide.^ It 
is a corrugated single layer of sp 2 hybridized hexagonal boron nitride with a lattice constant of 3.2 
nm. If the nitrogen atom sits close to an on top position of a substrate Rh atom, /z-BN is tightly 
bound and forms the so called "holes" or "pores" with a diameter of about 2 nm. The loosely bound 
regions form a connected network and are called wires. This structure imposes a template function 
which can be used for trapping single molecules at room temperature. 3 The trapping mechanism 
was found to be related to the presence of large lateral electric fields in the surface, where the rims 
of the holes carry lateral dipole rings, which produce a lower electrostatic potential in the holes.^ 
Also, it was shown that the /z-BN nanomesh is very robust and survives transport in air,^ immer- 
sion into liquids,^ and even electrolytic cycles.^ The impedance analysis in the electrochemical 
set up suggested that the /i-BN nanomesh imposes ordering in the electrolyte within the Helmholtz 
layer.^ This led us to the investigation of water on the nanomesh. It turns out that the hydrogen 
bonds dominate above the bonds due to the dipole rings. However, the rims of the holes prevent the 
ice clusters from further growth. It also becomes apparent that the wires act like a feeding/draining 
network that connects the ice clusters. 

Water on surfaces is an intriguing subject with many experimenta l 7 ^ 9 ! 1Q I 1 1 * 1 2 I and theoreti- 
ca || i3 | i4 | i5 | i6 | i7 | c ] ia n en g es j n i e t ter the /j-BN/Rh(lll) nanomesh is used as a template for the 
formation of ice clusters. It is found that at temperatures below 50 K water self-assembles in the 
nanomesh holes into two dimensional ice clusters with about 40 molecules, and remains as a low 
density phase on the wires. The structure of the clusters is elucidated, and in particular it is shown 
that the positions of the hydrogen atoms within the clusters may be inferred from tunneling barrier 
measurements with sub-nanometer resolution. 
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Experimental 

The experiments were performed in an ultrahigh- vacuum (UHV) chamber with a background pres- 
sure below 4xl0 -10 Torr using variable-temperature scanning tunneling microscopy (Omicron 
VT-STM). The sample preparation included several cycles of Ar + -bombardment, subsequent ex- 
posure to a few L (1 Langmuir=10~ 6 Torr-s) of O2 and annealing of the Rh(lll) sample. Then it 
was exposed to 40 L of borazine (HBNH)3 gas, while keeping the surface at 1070 K. This proce- 
dure yields a well-ordered large-scale single layer of hexagonal boron nitride on Rh(l 1 1) surface. 
Milli-Q water was used and purified by several freeze-and-pump cycles. Water was introduced 
into the UHV chamber via a leak valve in the pressure range of 10 10 -10 8 Torr through a nozzle 
pointing towards the sample. All pressures correspond to the uncorrected reading of the ion gauge 
of the UHV chamber. In order to minimize tip induced H2O motion on the surface, scanning pa- 
rameters were set to the order of 50 pA to 400 pA for tunneling currents and —50 mV to —400 
mV for tunneling voltages. Tunneling barrier spectroscopy was performed with lock-in technique 
by superimposing a modulation voltage of 8 mV at 2.0 kHz to the z-piezo. This modulates the 
tip-sample distance with an amplitude of ~0.07 nm. Except for those mentioned explicitly, STM 
images shown in this paper are obtained at low temperature 34 K. 



Results and discussion 



Figure 1 compares the /z-BN/Rh(lll) nanomesh with and without adsorbed water. The images 



with atomic resolution show that water forms a two-phase system, where the two phases coexist 



in the supercell with 3.2 nm lattice constant. Figure 1 1 is a large scale STM current image of the 
nanostructure (42x42 nm 2 ) across several terraces of the Rh(lll) substrate, obtained after dosing 
about 2 L water on the surface at 52 K. All holes are filled with ice clusters and no ordered water 
structures appear on the wires. This indicates that in the holes the interaction between water and 



the surface is strong enough to stabilize molecules. Compared to the empty mesh (Figure 1 :>), 



where the holes and the wires are visible as two structural elements of the nanomesh, the holes 
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are filled with small lobes and streaks appear along the scan direction after dosing 1 L of water 



< |Figure lb ). The topographic cross sections across a hole in |Figure Tp demonstrate that those water 
molecules map as protrusions with a height of about 15 pm, which is in agreement with the known 
mapping properties of water on surfaces. 1§ The lateral distance between the protrusions is about 
0.46 nm in the holes with a remaining depth of about 40 pm. As we show in the following, it is 
the signature of a bilayer of water — a lateral honeycomb structure with two sublattices A and 
B, where molecules in sublattice A form the contrast in the image. Furthermore, it is observed 
that no such ordered ice clusters appear on the wires, though they map rougher and noisy streaks 
reminiscent to mobile species are noticed. 



Figure 1 j shows that the height differences of the cuts around the wires are irregular, where 



the height variations are, compared to the empty nanomesh, larger for the water covered case. We 
speculate that, water molecules are moving fast along wires to form one low density gas phase 
which is e.g. similar to the 2D solid-gas structure obtained for SubPc molecules on Ag(lll) sur- 
face.^ These results indicate that the water self-assembly process in the boron nitride nanomesh, 
separates water into two distinct phases, one on the wires as fuzzy protrusions and one in the 
holes of the template as ordered nano-ice crystal. In the following we focus on the structure of the 
ordered ice clusters in the holes. 



The lattice constants and the azimuthal orientation can be seen in Figure 2 In Figure 2 1 and 



c, STM pictures of two different preparations are shown. The Fourier transforms (Figure 2 5 and 
d) reveal the two characteristic lattice constants of the nanomesh (3.2 nm) and the ice clusters 
(~0.46 nm). The lattice constant of ice is 7.0±0.2 times smaller than that of the nanomesh and 
corresponds within the error bar to that of the basal plane of hexagonal ice (0.45 nm). For the 
two shown preparations we find, different orientations rotated by 1 1±3° and -24±3° with respect 



to the direction given by two neighboring holes (Figure 2 3 and d). The observation of rigid ice 
clusters that are not oriented the same way with respect to the substrate indicate that the lock in 
energy of the nano-ice clusters is not as large as the bonding within the ice layer. The histogram in 
Figure 2b lists the number of molecular protrusions in the 32 holes of Figure 2 1 where in average 
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20±5 protrusions per hole were found. This confirms the above statement that water in sublattice 
A forms the contrast, since a hole with 2.2 nm diameter is expected to host twice as many i.e. about 
40 water molecules per bilayer. 



Figure 3 substantiates the two sublattice (bilayer) picture with a high resolution zoom into an 



ice cluster. The two sublattices A and B were resolved. We assign the highest protrusions A to 
positions of the topmost oxygen atoms, while the protrusions with intermediate height indicate 



oxygen atoms of sublattice B. From the cross section in Figure 3d it can be seen that sublattice 
B appears at the given tunneling conditions about 15 pm below sublattice A, where the direc- 
tion AB between nearest neighbor oxygen molecules encloses an angle of 93° with the surface 
normal. This confirms that at 34 K ice crystallizes in a bilayer structure in the holes of the h- 
BN/Rh(l 11) nanomesh. While it is intuitive to assign the two sublattices to two differently bound 
water molecules, the topographic images make no direct statement on the positions of the hydro- 
gen atoms, nor the orientation of the lonepairs on the oxygen atoms. In the following we show that 
the water dipoles and thus the positions of the hydrogen atoms affect the tunneling barrier. For this 
purpose models for the bilayer that obey the ice rules are considered. Essentially the ice rules say 
that every hydrogen atom forms a hydrogen bond with a neighboring oxygen atom. Every water 
molecule is modeled by a tetrahedron, with an oxygen atom in the center, two hydrogen atoms or 
protons in two of its vertices, and two oxygen lone pairs in the two remaining vertices of the tetra- 
hedron.^ For such a hexagonal bilayer the angle between AB and the surface normal is expected to 
be the tetrahedron angle of 109°. The deviation to the observed 93° is merely taken as an indication 
that the apparent heights in STM pictures must not correspond to the true topography, since it is 
also influenced by the local density of states and the tunneling barrier and that the tetrahedra might 
be distorted perpendicular to the surface. More importantly, the tetrahedra model does not make a 
statement on which vertices the protons sit. The basal hexagram plane is fully hydrogen bonded, 
i.e. half of the vertices are occupied with hydrogen atoms, and their nearest neighbor vertices are 
empty. It is worth to recall that the ice rules allow many different arrangements of the protons, 
i.e. allow for a certain proton disorder. The up (down) vertices either host a hydrogen (lone pair) 
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or vice versa. The issue on whether hydrogen atoms in sublattice A (B) point to the substrates 
or not, can be addressed by barrier height dl/dz microscopy ,120121] smce t h e positions of the ions 
leave fingerprints in the electrostatic landscape, as it was e.g. shown with image potential state and 
dl/dV spectroscopy for NaCl layers .^^Electrostatic potential variations locally affect the tun- 
neling barrier height. The tunneling barrier height is determined in measuring the z-dependence of 
the tunneling current. The tunneling current / corresponds exponentially on the distance z between 
the tip and the sample: / oc exp(— z/z ), where the tunneling length z depends on the tunneling 
barrier. For small tunneling voltages z is roughly proportional to 1/ V®i, where <!>/ is the local 
work function in the tunneling junction. In the barrier height microscopy mode the z-position of 
the tip is oscillated with a given frequency, while the tunneling current variations at this frequency 
are read out, and related via (dlnl/dz) 2 °« to the local work function. 20 Figure 4i shows a cor- 
responding dl/dz map of nano-ice//?-BN/Rh(l 1 1). Clearly, at the positions of the oxygen atoms of 
sublattice A and on the wires, a larger tunneling barrier is found, where there is some scatter in the 
barrier height above the individual up vertices. This is consistent with the assignment that the up 
vertices of sublattice A host oxygen lone pairs and correspondingly that the out of plane hydrogen 
atoms sit on the water molecules of sublattice B and point to the h-BN substrate. Density func- 
tional theory calculations of single water molecules on a single sheet of h-BN confirm this picture. 
They favor single water molecules to bind with the hydrogen to the negatively charged nitrogen 
atoms. 

For comparison we also show the dl/dz map of a bare nanomesh hole and the confining wires 
( |Figure 4p ). This result with atomic resolution on the /z-BN is in perfect agreement with the 



potential energy landscape of the clean nanomesh. 4 Both pictures in Figure 4i and b are recorded 
with the same tunneling parameters. The magnitude of the effect may be estimated from the 
measured potential energy difference 0.38 nm above the nitrogen atoms of 310 meV between 
holes and wires in the empty nanomesh.El 

In order to substantiate the dl/dz maps we built bilayer models of ice clusters obeying the 



above mentioned ice rules. A model of such an ice cluster is shown in Figure 4;. It shows 42 
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water molecules arranged in a honeycomb fashion, where in sublattice B the hydrogen atoms point 
to the substrate. There are several possibilities to arrange the hydrogen atoms. Here they are 
arranged in such a way that the cluster has no lateral dipolemoment. From density functional 
theory we take the Mullikan charges on the oxygen atoms of 0.4 electrons and —0.2 electrons on 
the hydrogen atoms (protons). With this we calculate the electrostatic potential above the cluster 
Figure 4b . Figure 4d shows this potential 0.2 nm above the plane between sublattice A and B. 



in 



Clearly, the potential that is reminiscent to a local work function is largest on top of the oxygen 
atoms in sublattice A, i.e. on top of the lone pairs that are not saturated with hydrogen atoms. For 
the chosen parameters we find a corrugation of the electrostatic potential energy landscape of up 
to 0.8 eV within the cluster. It has to be noted that this energy corrugation is distance dependent, 
and reproduces the order of magnitude of the effect. The calculations of the tunneling barrier 
also indicate that the dl/dz measurements feature the electrostatic potential closer to the surface 
than adsorbed Xe (0.38 nm), 4 since at this distance the molecular resolution is lost. All other 
configurations we modeled, like e.g. a structure with hydrogen atoms pointing to the vacuum, do 
not produce this kind of agreement. It is interesting to point out that the variations of the tunneling 
barrier heights on the different water molecules in sublattice A indicate that the method is sensitive 
to the local arrangement of the dipoles of the water molecules and that this bears the potential to 
address the positions of the protons. Also the findings indicate that different proton configurations 
near the edge of the nano-ice clusters cause electric fields that are important for the understanding 
of the growth limitation of the clusters and their reaction against other host molecules. 

In conclusion we report on the self-assembly of nano-ice clusters with about 40 molecules and a 
dilute "gas phase" on the /z-BN/Rh(l 1 1) nanomesh template. The nanomesh holes host the clusters 
and limit their size. With the help of barrier height dl/dz spectroscopy the structure of the nano-ice 
could be refined to that of a hydrogen bonded honeycomb bilayer, where the out of plane hydrogen 
atoms bond toward the surface. Also the positions of the protons in the layer may be inferred from 
the orientations of the water dipole fields. This new structure of ice gives unprecedented insight 
into the self assembly process of water and indicates it to be a good candidate to study proton 
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Figure 1: Nano-ice on the h-BN/Rh(lll) surface, (a) Large-scale STM current image of water 
molecules adsorbed on the surface at T=52 K. All holes are occupied by distinct ice clusters (42 
nmx42 nm; Y s = -0.4 V; L= 400 pA). (b) Atomically-resolved STM topography images of the bare 
surface and after dosing water (c), both areas are 9 nmx9 nm; and scanning parameters are Y s = 
-0.05 V; L= 100 pA. Height profiles across (d) and around (e) the holes in (b) and (c) (blue lines 
and red lines, respectively). The height level in (d) corresponds to the average of (c) and that in 
(e) to its average. 
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Figure 2: Two different preparations of nano-ice clusters on the surface, (a) and (c) STM topogra- 
phy images of nano-ice on h-BN/Rh(l 1 1) nanomesh (18 nmx 18 nm). (b) and (d) The correspond- 
ing Fourier transforms of (a) and (c), respectively. The contours of a Gaussian anchor the exact 
position for the periodicity of water molecules in reciprocal space. The inset angles corresponding 
to the orientations between the direction of nano-ice and its neighboring holes, (e) The histogram 
of the numbers of top layer H2O molecules per hole for the area in (a). 
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Figure 3: Atomic-scale structure of the nano-ice cluster revealing the honeycomb bilayer with two 
sublattices A and B. (a) High-resolution STM image with the positions of the oxygen atoms (red 
balls A and B) are embedded (0.8 nmx 1.2 nm; V J = +0.3 V; l t = 100 pA). (b) Height profile across 
the water molecules in (a) as indicated by the yellow line. The two different water molecules are 
schematically indicated. The angle between AB and the surface normal is 93°. 
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Figure 4: (a) The dl/dz map of nano-ice cluster in the hole after dosing 1L H2O molecules on the 
surface. By measuring the response of the tunneling current to the variation of the tunneling gap 
distance, dl/dz image can be obtained simultaneously with a topographic STM image. For the 
dl/dz map, the bright colors correspond to large local tunneling barriers, (b) The dl/dz map was 
measured on the bare surface (Both images are 6 nmx6 nm; scanning parameters are Y s = -0.05 
V; l t - 100 pA). (c) Side and top views of the structure model for nano-ice cluster in a single hole, 
(d) The Corresponding surface electrostatic potential 0.2 nm above the plane between sublattice A 
and B in (c) (3 nmx3 nm). 
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